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In order to probe the effect of charge fluctuations on triplet pairing, we study the pairing sym- 
metry in the one-band Hubbard model having the off-site Coulomb repulsion (V) on top of the 
on-site repulsion as a model for the 7 band of SkRuCU, a strong candidate for triplet pairing su- 
perconductor. The result, obtained with the dynamical cluster approximation combined with the 
quantum Monte-Carlo method, and confirmed from the fluctuation exchange approximation, shows 
that while d I .2_ !/ 2-pairing dominates over p in the absence of V, introduction of V makes p x +y and 
d xy dominant. The gap function for the chiral p x + y -\-ip x - y has nodes that are consistent with the 
recent measurement of specific heat in rotated magnetic fields in the ruthenate. This suggests that 
the off-site repulsion may play an essential role in triplet superconductivity in this material. 

PACS numbers: 74.70.Pq, 74.20.Rp, 71.10.Fd 



Introduction — There has been an increasing fasci- 
nation with spin-triplet pairing in the condensed mat- 
ter physics, with a history dating back to the discovery 
of superfiuid 3 He. While triplet superconductivity has 
been found in some heavy fermion compounds and or- 
ganic metals, a ruthenate, Sr2Ru04, is prototypical in 
that its structure is similar to a cuprate La2Cu04, but 
the replacement of Cu with Ru makes the relevant d or- 
bitals different. So the discovery of superconductivity^ 
in the ruthenate has stimulated an enormous amount 
of studies^]- Experimentally, it has been established 
that Sr2RuC>4 is a chiral (time-reversal broken) triplet 
superconductor with the spins lying in the RuC>2 plane, 
where we take the spin quantization axis (z) along the 
crystalline c-axis, as indicated from Knight shift 3| and 
/xSR0] experiments. 

On the other hand, the symmetry of the gap function, 
which is of prime importance in identifying the pairing 
mechanism, has yet to be established. Although the exis- 
tence of line-nodes in the gap function is suggested from 
power-law behaviors in specific heat and NMR l/TiT@, 
the position of nodes remains controversial. While the 
magnetothermal conductivity shows a weak anisotropy 
in the a6-plane for T four-fold oscilla- 

tions [indicative of nodes in the gap function around 
k =(±7r,0), (0,±7r)] have been detected in a recent mea- 
surement of specific heat in rotating magnetic fields by 
Deguchi et aZ.pfl. 

One complication is that, unlike the high-Tc cuprates 
which have a square lattice of d x 2_ y 2 orbitals, the ruthen- 
ate has three, cylindrical Fermi surfaces (labelled as a, 
(3, and 7), where a and j3 derive from one-dimensional 
arrays of Ru d xz and d yz orbits, while 7 derives from a 
square lattice of d xy orbits. This has been established 
experimentally from a de Haas measurement and an 
angle-resolved photoemission spectroscopy, and also 
confirmed from a first-principles band calculation|ld|. 



Deguchi et al. have concluded that their result sug- 
gests that the active band for superconductivity is the 
7 band.0 

Although various pairing mechanisms have been pro- 
posed for Sr2Ru04, they are not straightforward, since 
it is much more difficult to explain triplet pairing than 
singlet pairing solely from spin fluctuations because of a 
smaller pairing interaction in the triplet channel. Sev- 
eral authors 0, 0, 0, Q have focused on the effect 
of nesting in the quasi-one-dimensional a-(3 Fermi sur- 
faces to show that anisotropic antiferromagnetic spin 
fluctuations, observed in NMR experiment [15J, or or- 
bital fluctuations can favor the triplet pairing. On 
the other hand, assuming that the 7 band is the ac- 
tive band for superconductivity, a one-loop renormal- 
ization group analysis for the one-band Hubbard model 
was performed 1 16| . where p-pairing is concluded to be 
dominant unless the Fermi level is far from the van 
Hove singularity. A third-order perturbation calculation 
for the single-band Hubbard model has also been per- 
formedei q/.|l7t Il8| . where triplet pairing is shown to 
dominate over singlet pairing for intermediate band fill- 
ing. This was recently extended to the three-band Hub- 
bard model. [l^. fig|| for which the perturbation calcula- 
tion shows that triplet pairing remains dominant, resid- 
ing on the 7 band. Their result that the gap function 
has nodes along [100] for the 7 band, while they lie along 
[110] for a- (3, is consistent with the anisotropy in the 
magnetothermal conductivity 0, However, we believe 
that the validity of the results obtained with perturba- 
tion expansions, truncated at the third or fourth order, 
has to be checked by non-perturbative methods. 

Indeed, according to the fluctuation exchange (FLEX) 
study [i(J or a phenomenology [2]] for the spin fluctua- 
tion, triplet pairing is rather weak in general. More 
specifically, Kuroki et al. have shown recently that the 
singlet pairing dominates over the triplet for the one- 
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band Hubbard model for the 7 band, where they cal- 
culated the pairing interaction vertex with the quantum 
Monte Carlo (QMC) method to show that d-pairing in- 
ter action is stronger than those in triplet channels. |22| 

Given this background, the motivation of the present 
study is two- fold. First, we want to clarify whether the 
triplet superconductivity can be dominant within the 
one-band Hubbard model in a non-perturbative method, 
for which we have employed the dynamical cluster ap- 
proximation (DCA) combined with QMC formulated by 
Jarrell et of. [23| . We shall show that d x 2_ y 2-w&ve pairing 
dominates over p-w&ve pairings suggested by third-order 
perturbation studies. Here DCA+QMC is employed for 
the first time[27| to investigate superconductivity in the 
extended Hubbard Hamiltonian that models the 7 band 
in the ruthenate. 

Second, and more importantly, we propose that the 
off-site Coulomb repulsion can favor triplet supercon- 
ductivity in general and in the 7 band in Sr2Ru04 in 
particular |2J|. The underlying physics is that, while we 
usually evoke spin fluctuations (caused by the on-site re- 
pulsion) in considering an electron mechanism of super- 
conductivity, charge fluctuations, which tend to be en- 
hanced by the off-site repulsion, should favor the pairing 
in the triplet channel, as seen in the expression for the 
fluctuation-mediated interactions 0, H^}- We have in 
fact shown with a FLEX calculation that, even for the 
simple square lattice, transitions between different pair- 
ing symmetries can arise due to the coexistence of spin 
and charge fluctuations [2(|. We shall show that the off- 
site repulsion makes p x +y and d xy dominant while p x and 
d x 2_ y 2 are suppressed. The result is also confirmed by 
a FLEX calculation. The chiral p x + y +ip x - y gap func- 
tion then has nodes around (±7r,0) and (0,±7r), which 
is consistent with the specific heat in rotating magnetic 
fields 0. 

Formulation — The extended Hubbard model is given 

as 

nn nnn 

n = -tJ2 kW + h.c.)-t' J2 kW + h.c.) 

nn 

+U n^riji + V n ia n. j<yl 

in the standard notation, where nn (nnn) denotes 
nearest-neighbor (next-nn) sites. In DCA, the original 
reciprocal space with TV points are divided into N c cells, 
for which a coarse-graining is done. The Hubbard model 
is then mapped to a self-consistently embedded cluster of 
Anderson impurities (rather than a single impurity con- 
sidered in the dynamical mean- field approximation 28J), 
so that DCA incorporates nonlocal spatial fluctuations. 
We have solved the cluster problem generated by the 
DCA using the QMC with the algorithm proposed by 
Hirsch and Fye 29] . We choose a cluster size N c = 4 x 4 
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FIG. 1: (a) DCA+QMC result for the pairing susceptibil- 
ity versus temperature for various symmetries of pairing (de- 
picted in real space in inset) in the on-site Hubbard model 
(V — 0) with the on-site interaction U = 1.5, the band filling 
n = 1.33, and t' /t — 0.4. Error bars are smaller than the size 
of each symbol. 



throughout the study [30(. 

The effect of the off-site repulsion V can be incorpo- 
rated into the DCA calculation as follows. The off-site 
term is expressed in fc-space as 
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with Vq = (V/N)[cos(q x ) + cos(q y )], which has to be 
coarse-grained. One might think that the coarse-graining 
may introduce interactions extending beyond nearest 
neighbors for the impurity cluster model. However, 
we can note a relation, Vq = (N c /N) J2q ^Q+q = 
(V /7V c )(cos Q x +cos Q y ), where the summation for q runs 
over the momenta in the coarse-graining cell centered at 
Q, and V = sm(ir/N c )/(ir/N c )V. Thus the form of the 
interaction does not change, so we have only to consider 
the off-site repulsion V in the QMC for the impurity clus- 
ter model. 

Results — In the DCA+QMC calculation, we take the 
half of the band width as the unit of energy as customary 
done. We take the on-site Coulomb interaction U = 1.5, 
the band filling n = 1.33, and t'/t = 0.4. These val- 
ues are chosen to roughly represent the 7 band of the 
ruthenate. The pairing symmetries considered are s *~ 1, 
p x - v / 2sin(fc a; ), 

Px±y ~ \/2sin(k x i k y ), 31] d x 2_ y 2 ~ 
cos(k x ) — cos(k y ), and d xy ~ cos(k x + k y ) — cos(k x — k v ), 
as depicted in real space in the inset of Fig^ We first 
show the result before V is turned on. Figure^shows the 
pairing susceptibility for the on-site Hubbard model as a 
function of temperature, where L = 64 Trotter-Suzuki 
decomposition number is taken. We can see that d x 2_ y 2 
is the dominant pairing, while p x + y is weaker in the re- 
gion T > 0.031 studied here. The result is consistent with 
Ref. where d x 2_ y 2 is shown to dominate over triplet 
pairings by a QMC for a finite (N = 14 x 14) Hubbard 
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FIG. 2: The pairing susceptibility in the extended Hubbard 
model as a function of the off-site repulsion V for U = 1.5, 
n = 1.33, and t /t = 0.4 at T = 0.125. Error bars are smaller 
than the size of each symbol. 

model with U = 0.5. These results, obtained with non- 
perturbative methods, suggest that the on-site Hubbard 
model is insufficient to describe the superconductivity in 

Sl'2Ru04. 

Now we move on to the extended Hubbard model in 
Fig|21 where the DCA+QMC result for the pairing sus- 
ceptibility is plotted as a function of the off-site interac- 
tion V for T = 0.125 with the Trotter-Suzuki decompo- 
sition number L — 24. We can clearly see a qualitative 
tendency that p x + y and d xy become dominant with the 
introduction of V, while p x and d x 2_ y v are suppressed. 
Physically, this should be because the nearest-neighbor 
V suppresses pairs formed across nearest-neighbor sites, 
while the pairs such as p x + y and d xy that are formed 
across more distance sites (see the inset of Fig^) are 
less affected. We cannot take a lower T or larger L 
at present, since in the QMC algorithm for V^nia-riia-i 
eight Hubbard-Stratonovich auxiliary fields (for er, a' =\ 
, J, for x, y directions) are needed on top of the one for the 
[/-term, so the study of the subtle competition between 
p x + y and d xy at lower temperatures is a future study [32|. 

Still, it is interesting to discuss what should result from 
the p-wave pairing. Below T c the p x + y -wetve pairing is 
expected to take a chiral form, 

p x + y + i>Px-y ~ sin(k x + k y ) +i sin(k x - ky), 

since the mixing should increase the gap energy (| A|). We 
can then note, as depicted in Fig|3 that the gap function 
A for the p x + y + ip x - y has nodes (or minima of |A| on 
the Fermi surface) at around (±7r,0) and (0,±tt) (which 
happen to be similar to those for the chiral p x + ip y ~ 
siri(k x ) + i sin(fc a ); FigOJb)). So the chiral p x + y + ip x - y 
is consistent with the specific heat measurement Q, and 
should be a candidate for the triplet superconductivity 
in S^RuO^. 

We finally examine whether the FLEX, a renormal- 
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FIG. 3: (a) The position of nodes of p x + y ~ sin(k x + k y ) (solid 
lines), p x - y ~ sm{k x - k y ) (dotted), and p x+y + ip x - y (gray 
dots). The bold solid curve represents the 7 Fermi surface of 
Sr2Ru04. (b) A plot similar to (a) for p x + ip y ~ sin(k x ) + 
i sin^y). 
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FIG. 4: Contour plot of the gap function for the singlet(a) and 
triplet (b) pairs obtained with FLEX+Eliashberg's equation 
for the extended Hubbard model for U = 2.5, V = 0.65, 
n = 1.33, t'/t = 0.5, and T = 5 x 10~ 3 . White dashed lines 
denote the nodes. For (b), another one, rotated by 90 degrees, 
is degenerate. 

ized perturbation, would reproduce the above result. The 
FLEX was first formulated by Bickers et al.^ij^ for the 
Hubbard model, and further applied to the extended 
Hubbard model by Esirgen et al. [34| . We first obtain 
the renormalized Green's function, G, taking bubble and 
ladder diagrams as the self energy. We then calculate the 
pairing interaction mediated by spin and charge fluctua- 
tions, and plug that in Eliashberg's equation, 

WW = -§E rPP ( fc - fc ')G(fc')G(-fcW), (1) 
fe' 

where <fi is the gap function, k = (k,uj n ) with ui n being 
Matsubara frequency, and r pp the pairing interaction be- 
tween the pairs with (fc, —k) and (fe , — fe ). 

The maximum eigenvalue A becomes unity when T 
becomes the transition temperature of the dominant 
pairing. We take N = 32 x 32 sites at temperature 
T = 5 x 10~ 3 and -(2N C - l)vrT < u>„ < (2N C - 1)ttT 
with N c = 1024. When the off-site interaction V is 
introduced all the vertices and susceptibilities become 
[Z + 1) X (Z + 1) matrices for the lattice coordination 
number Z(= 4 for the square lattice). 

We have performed the FLEX calculation for t'/t = 
0.4,0.5, n = 1.33, U = 1.5 - 2.5, V = - U/4. While 
the singlet pairing is dominant in most of the parameter 
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region, we have found that the maximum eigenvalue of 
Eliashberg's equation becomes greater for triplet pairing 
than that of the singlet for 0.625 < V < 0.675 for U = 
2.5|33. Figure0]depicts the gap functions for triplet and 
singlet cases. The symmetry of the gap function, which 
is Px+y- and d xy -like, res pec tively, does agree with the 
above DCA+QMC result. 

In summary, we have found using DCA+QMC and 
FLEX methods that, while p-wave pairing is weaker than 
d x 2_ y 2 in the on-site Hubbard model, introduction of the 
off-site repulsion V suppresses p x and d x 2_ y 2, making 
p x + y and d xy more favorable. Although the effect of a 
and (3 bands is beyond the scope of the present study, the 
qualitative tendency of the present work that the off-site 
repulsion acts to preferentially suppress nearest neigh- 
bor pairs should hold for multibands as well as for single 
bands. The position of the nodes in the gap function for 
the chiral p = p x + y + ip x ~ y is consistent with a specific 
heat measurement for Sr2RuO4 0, so this is a promising 
candidate for the triplet superconductivity in Sr2Ru04, 
and the off-site repulsion may play an essential role there. 

This work was supported by a Grant-in-aid for scien- 
tific research and Special coordination funds from the 
ministry of education of Japan. Numerical calculations 
were performed on SR8000 in ISSP, University of Tokyo. 
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